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In this paper we build a holographic model of paramagnetism/antiferromagnetism phase transi- 
tion, which is realized by introducing two real antisymmetric tensor fields coupling to the background 
gauge field strength and interacting with each other in a dyonic black brane background. In the case 
without external magnetic field and in low temperatures, the magnetic moments condense sponta- 
neously in antiparallel manner with the same magnitude and the time reversal symmetry is also 
broken spontaneously (if boundary spatial dimension is more than 2, spatial rotational symmetry 
is broken spontaneously as well), which leads to an antiferromagnetic phase. In the case with weak 
external magnetic field, the magnetic susceptibility density has a peak at the critical temperature 
and satisfies the Curie-Weiss law in the paramagnetic phase of antiferromagnetism. In the strong 
external magnetic field case, there is a critical magnetic field Be in antiferromagnetic phase: when 
magnetic field reaches Be, the system will return into the paramagnetic phase by a second order 


phase transition. 


I. INTRODUCTION 


The anti-de Sitter/conformal field theory (AdS/CFT) 
correspondence relates a weakly coupled gravity theory 
in anti-de Sitter (AdS) space to a strongly coupled con- 
formal field theory (CFT) living on the AdS boundary [1- 
4]. Due to the existence of scaling symmetry near critical 
point, the AdS/CFT correspondence provides a powerful 
tool to study critical phenomena in some strongly cou- 
pled condensed matter systems. This duality also gives 
us a way to understand gravity and condensed matter 
physics from other side. Over the past years there have 
been a lot of studies of the applications of the AdS/CFT 
duality in condensed matter physics (for reviews, see 
Refs. [5-9]). Recently, some efforts have been made to 
generalize the correspondence to systems with less sym- 
metries (see Refs. [10-16], for example) and to the far 
from thermal equilibrium problems (see Refs. [17-21], for 
example). 

In a previous paper [22], the present authors proposed 
a new example of the application of the AdS/CFT cor- 
respondence, by realizing the holographic description of 
the paramagnetism/ferromagnetism phase transition in 
a dyonic Reissner-Nordstr6m-AdS black brane. In that 
model, the magnetic moment is realized by condensa- 
tion of a real antisymmetric tensor field which couples 
to the background gauge field strength in the bulk. In 
the case without external magnetic filed, the time rever- 
sal symmetry is broken spontaneously and spontaneous 
magnetization happens in low temperatures. The critical 
exponents are in agreement with the ones from mean field 
theory. In the case of nonzero magnetic field, the model 
realizes the hysteresis loop of single magnetic domain and 
the magnetic susceptibility satisfies the Curie-Weiss law. 

Except for the paramagnetism and ferromagnetism, 
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antiferromagnetism is another important magnetic prop- 
erty of material. So a natural question is whether we 
can extend the model in Ref. [22] to realize antiferro- 
magnetism. To do that, let us first briefly review some 
characteristic properties and some relevant theory of an- 
tiferromagnetism [23, 24]. 

Though an antiferromagnetic material does not show 
any macroscopic magnetic moment when external mag- 
netic field is absent, it is still a kind of magnetic ordered 
material when temperature is below the Néel tempera- 
ture Ty. The conventional picture, due to L. Néel, rep- 
resents a macroscopic antiferromagnetism as consisting 
of two sublattices 1, call them A and B such that spins 
on one sublattice point opposite to those of the other 


sublattice 7, i.e., Ma = —Mp. The order parameter 
is the staggered magnetization, defined as the difference 
between the two magnetic moments associated with the 
two sublattices: 


Mt =Ma- Mp. (1) 


Such a complicated ordering is not straightforward to 
identify by thermodynamic measurements, and in early 
days even the existence of antiferromagnetism was hotly 
debated. However, the technology of neutron diffraction 
gives us a powerful tool to directly measure the micro- 
scopic magnetic moment distribution of the lattices. The 


1 There are some antiferromagnetic materials which have more 
than two sublattices. In those materials, the spins in the sublat- 
tices are arranged with zero total spin. These kinds of materials 
can be regarded as the extension of the two sublattices. In this 
paper, as a toy model, we will considere the case with two sub- 
lattices only. 

2 There are some antiferromagnetic materials, in which the mag- 
netic moments are a little canted with each other near the tran- 
sition temperature, such as œ-Fe203 in the temperature of 950- 
260K. This phenomenon originates from some scattering and 
perturbation effects of spin-orbit coupling [25]. However, when 
temperature approaches to zero, the magnetic moments become 
antiparallel. In this paper, we will not consider these cases. 


experiments confirmed this picture in the antiferromag- 
netic insulator. Figure I depicts the spin alignments for 
two manganese compounds. (Only the spins of the man- 
ganese ions contribute to the antiferromagnetic behav- 
ior.) Figure I (a) shows that the ions in a given {111} 
plane possess parallel spin alignment, whereas ions in the 
adjacent plane have antiparallel spins with respect to the 
first plane. In Figure I (b), the spins of manganese ions 
in vertexes are all parallel alignment, whereas ions in the 
body-centered are arranged in antiparallel manner with 
vertexes. Thus, both in these two materials, the magnetic 
moments of the solid cancel each other and the material 
as a whole has no net magnetic moment. 


(a) 


(b) 


FIG. 1. Schematic representation of spin alignments for an- 
tiferromagnetics at 0 K. (a) Display of magnetic structure in 
MnO. The arrows stands for the directions of spin and circles 
represent the Mn ions [26]. The oxygen ions (are not shown 
here) do not contribute to the antiferromagnetic behavior. 
MnO has a NaCl structure. (b) Three-dimensional represen- 
tation of the spin alignment of manganese ions in MnF»2 [27]. 
The smaller spheres stand for Mn ions and larger ones stand 
for fluoride ions. This figure demonstrates the interpenetra- 
tion of two manganese sub-lattices, Mn; and Mnog, having 
antiparallel aligned moments. 


One can see that, below the Néel temperature, the 
magnetic moment in antiferromagnetic material is still 
ordered rather than the case in paramagnetism where 
the magnetic moments are arranged randomly because 
of thermal fluctuations. However, when the tempera- 
ture increases beyond Ty, the ordered magnetic moment 
structure is broken and the magnetic moment begins to 
distribute randomly, which corresponds to the paramag- 
netic phase of the antiferromagnetic material. 

There is another significant difference between the an- 
tiferromagnetism and paramagnetism which can be seen 
from magnetic susceptibility. The behavior of magnetic 
susceptibility of antiferromagnetism is shown schemati- 


cally in figure 2. In the high temperature region T > 
Ty, the antiferromagnetic material is in a paramag- 
netic phase whose magnetic moments are randomly dis- 
tributed. So the susceptibility performs as a paramag- 
netism and obeys the Curie-Weiss law [28]: 
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where C and @ are two constants. Note that in the 
paramagnetic material and the paramagnetic phase of 
ferromagnetic material, the magnetic susceptibility also 
obeys the Curie-Weiss law, but the constant @ in (2) is 
zero and negative, respectively. For the antiferromag- 
netic material, in the region of low temperature T < Ty, 
the magnetic moments condense and arrange into two 
opposite directions. In an external magnetic field, a kind 
of ferrimagnetic behavior may be displayed in the anti- 
ferromagnetic phase, with the absolute value of one of 
the sublattice magnetizations differing from that of the 
other sublattice, resulting in a nonzero net magnetiza- 
tion. For a weak external field, there are two cases ac- 
cording to the directions of external magnetic field and 
spontaneous magnetization. When the external field is 
perpendicular to spontaneously induced magnetic mo- 
ment, the magnetic susceptibility is almost independent 
of temperature, while in the parallel case, the suscep- 
tibility decreases when temperature is lowered, because 
the force making magnetic moments antiparallel becomes 
stronger and stronger, which leads to that the material 
is harder and harder to be magnetized by external field. 
However, if the external field is enough strong, it will de- 
stroy the antiferromagnetic interaction of the two mag- 
netic moment. Then the system will perform a param- 
agnetism, i.e., both two sublattices have same magnetic 
moments paralleling the external magnetic field as if they 
have no interaction with each other (see the right plot of 
figure 2). 

In general, magnetism of materials originates from the 
exchange interaction of electrons. There are two meth- 
ods to deal with the magnetic ordering in materials. One 
is based on the local magnetic moment and the basic 
model is Heisenberg Hamiltonian. The other is based 
on the energy band theory with the effects of spin and 
the basic model is Hubbard Hamiltonian. The former is 
suitable to describe magnetic properties in insulator and 
some semiconductor, while the latter is suitable for con- 
ductor. In both models, the exchange interaction plays 
a crucial role. Various microscopic exchange interactions 
between the magnetic moments or spins may lead to anti- 
ferromagnetic structures. Generally speaking, a negative 
exchange interaction is the main reason for antiferromag- 
netism. In the simplest case, one may consider an Ising 
model on an bipartite lattice, e.g., the simple cubic lat- 
tice, with couplings between spins at nearest neighbor 
sites. Depending on the sign of the interaction, ferromag- 
netic or antiferromagnetic order will appear. Geometri- 
cal frustration or competing between ferro- and antiferro- 
magnetic interactions may lead to different and, perhaps, 
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FIG. 2. (Top panel) The magnetic susceptibility of antifer- 
romagnetic material MnF»2 [24]. xı stands for the case with 
external magnetic field perpendicular to the direction of spon- 
taneous magnetization and x for the case with external mag- 
netic field parallel to the direction of spontaneous magnetiza- 
tion. (Bottom panel) The influence of external magnetic field 
on the magnetic moments of two sublattices when T < Ty. 


more complicated magnetic structures [28]. These theo- 
ries give us a very deep understanding in the weak cor- 
relation system. However, in recent years, many antifer- 
romagnetic or ferromagnetic materials show that strong 
correlation plays important role and exhibits some novel 
properties [29-31]. So for these materials, the AdS/CFT 
correspondence might provide a useful method to under- 
stand these properties. 

The typical paramagnetism/antiferromagnetism phase 
transition is a second ordered phase transition. So this 
phase transition should be associated with some symme- 
try breaking spontaneously. Normally the appearance 
of antiferromagnetism (or ferromagnetism) is associated 
with the breaking of a continuous symmetry-spatial rota- 
tions and a Zz symmetry-time reversal simultaneously in 
three dimensions®. Though the former is more obvious, 
the latter is a characteristic feature. Since the magnetic 
field and spin (or angular moment) are vectors (in fact 
they are axial vectors) only in three spatial dimensions. 
As a theoretical consideration, we can also build an anti- 
ferromagnetism (or ferromagnetism) models in two spa- 
tial dimensions*, where the magnetic field and spin (or 
angular moment) are pseudo scalars. In these models, 
the magnetic ordered phase only breaks the time rever- 
sal symmetry spontaneously. In a strict sense, if a field or 


3 The magnetic moment is proportional to the spins and orbital 
angular moment (in some case, the latter can be neglected), so it 
will obtain a minus sign under the time reversal transformation. 

4 In two spatial dimension, isotropic Heisenberg model with finite- 
range exchange interaction can be neither ferromagnetic nor an- 
tiferromagnetic [32]. However, it does’t exclude anisotropic cases 
or the models beyond Heisenberg model. 


operator O is invariant under the time reversal transfor- 
mation, there is no interaction such as BO with respect 
to the external magnetic field B. 


Note that a holographical antiferromagnetic phase is 
proposed by breaking a global SU(2) symmetry repre- 
senting spin into a U(1) subgroup in a charged black hole 
background in Ref. [33]. The symmetry breaking is trig- 
gered by condensation of a triplet scalar field charged 
under the SU(2) gauge field. Though SO(3) symme- 
try can be transformed into SU(2) symmetry and this 
model leads to the spatial rotational symmetry breaking 
spontaneously, the time reversal symmetry is not broken 
spontaneously in the magnetic ordered phase. On the 
other hand, the attention in [33] is focused on the sym- 
metry breaking of paramagnetism/antiferromagnetism 
phase transition and the dispersion relation of the mag- 
netic excitations. The structure of magnetic moment 
in material and the response to external magnetic field, 
which are the characteristic properties of antiferromag- 
netism, are not manifest in that model. So from the point 
of view of symmetry breaking and experimental phe- 
nomenological analysis, Ref. [33] seemingly does not give 
a complete model for paramagnetic/antiferromagentic 
phase transition. 


Therefore a holographic model which is said to real- 
ize the phase transition of paramagnetism and antiferro- 
magnetism should at least give three features: (1) The 
antiparallel magnetic moments structure as T < Ty. (2) 
The susceptibility behavior shown in figure 2. (3) Break- 
ing of time reversal symmetry, and of spatial rotational 
symmetry if spatial dimension is more than 2. These 
are just what we want to realize in this paper. We will 
generalize the model in Ref. [22] to describe the para- 
magnetism/antiferromagnetism phase transition in a dy- 
onic black brane background by introducing two real an- 
tisymmetric tensor fields which couple to the background 
gauge field strength and interact with each other in the 
bulk. In the case without external magnetic field and 
in low temperatures, the magnetic moments of these two 
tensors condense spontaneously in the antiparallel man- 
ner with the same magnitude, which leads to an anti- 
ferromagnetic phase. Since these two magnetic moment 
both obtain minus signs under the time reversal trans- 
formation, the time reversal symmetry is broken sponta- 
neously. If the boundary spatial dimension is 3, magnetic 
moment is an axial vector and the spatial rotational sym- 
metry is broken spontaneously as well. In the case with 
weak external magnetic field, the magnetic susceptibility 
has a peak at the critical temperature and satisfies the 
Curie-Weiss law in the paramagnetic phase of antiferro- 
magnetic materials. For strong external magnetic field, 
there is a critical magnetic field Be. If B > Be the two 
magnetic moments in sublattices will have the same value 
paralleling to external magnetic field. 


The paper is organized as follows. In section II, we 
will first introduce the model and derive its equations of 
motion (EoMs). Then we will compute the on shell Eu- 
clidean action of the model, which is equivalent to give 


the grand potential of dual boundary theory. The numer- 
ical calculations on the magnetic moment structure and 
the susceptibility in different phases will be given in sec- 
tion III. Because the calculations in the 3+1 bulk dimen- 
sional and 4+1 bulk dimensional cases are very similar, 
we will consider the 3+1 bulk dimensional case for sim- 
plification ° and give some comments for the 4+1 dimen- 
sional case. We will point out some differences and give 
some discussions in section IV for the 4+1 dimensional 
case. The summary and some discussions are included in 
section V. A description in terms of mean field theory for 
the antiferromagnetism/paramagnetism phase transition 
is given in appendix A. 


Il. HOLOGRAPHIC MODEL 


A. Model and EoMs 


Generalizing the model proposed in [22], let us consider 
the following action: 


1 
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L is the radius of AdS space, 2k? = 167G with G the 
Newton constant. In this model, k, m? and J are all 
model parameters with J < 0. \? characterizes the back 
reaction of the two polarization fields M, A) with a = 1,2 
to the background geometry, and L12 describes the inter- 
action between two polarization fields. When k = 0, the 
two polarization fields decouple. In that case, the two 
polarization fields condense independently at the same 
critical temperature with the same magnitude. However, 
when k Æ 0, because of the interaction between them, re- 
sulting magnetic moments are dependent. In some suit- 
able range of the model parameters, two condensed mag- 
netic moments may appear in an antiparallel manner. 
Since two rank-two fields involve in this model, some 
remarks are in order. It is a well-known fact that it is 
quite difficult to write down a self-consistent high spin 
field theory in a flat spacetime (see [34], for example). 


5 In AdS/CFT duality, 3+1 gravity is dual to a 2+1 boundary 
theory, which can be regarded as an effective description for some 
antiferromagnetic film in real materials. 


Usually high spin field theory suffers from ghosts and 
causality violation. But such problems do not exist in 
this model and the ferromagnetic model in Ref. [22]. 
The reason is that our rank-two fields are two antisym- 
metric fields, which in fact are two massive 2-form fields 
rather than two spin-2 fields. The latter requires the 
field to be a symmetric and traceless one. As a result, 
the spin degree in general dose not agree with the rank 
of the field and depends on the dynamic of the field [35]. 
For example, the 2-form Kalb-Ramond field is spin O for 
massless case [36] and spin 1 for massive case [37]. p-form 
field in string/M-theory is well defined, as the source of 
the (p — 1)-brane. We can obtain a well defined theory 
for any form field localized in D-branes [38], which is en- 
tirely different from the case involving symmetric high 
rank tensor field. The spin of standard p-form field in 
type IIB supergravity theory can be found in the page 76 
of Ref. [35], where one can see that p-form field does not 
suffer from the troubles of high spin field. The mass of 
p-form field can be generated through Higgs mechanism, 
Stueckelberg mechanism, or topological mass generation 
method. As a result, the model here and the one pro- 
posed in Ref. [22] do not involve the higher spin fields 
and we do not worry the troubles such as ghost, causality 
violation and so on. 
The equations of motion for polarization fields read 
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In the probe limit, we can neglect the back reaction of 
the two polarization fields on the background geometry. 
The background we will consider is a dyonic Reissner- 
Nordström-AdS black brane solution of the Einstein- 
Maxwell theory with a negative cosmological constant, 
and the metric reads [39] 
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Here both the black brane horizon and AdS radius have 
been set to be unity. The temperature of the black brane 
is 

1 


T= 
At 


(3— p’ — B?) (7) 
For the solution (6), the corresponding gauge potential 
is A, = (1 — 1/r)dt + Bady. Here u is the chemical 
potential and B can be reviewed as the external magnetic 
field of dual boundary field theory. 

In the background (6), let us consider the dynamics of 
the two polarization fields. Following [22], we take the 


ansatz for the polarization fields M, (a) as 


M = —p™® (r)dt A dr + p® (r)dz A dy, 


a=1,2. (8) 


For the sake of convenience for later discussions, we in- 
troduce two variables 


1 1 
a= slo +p), B=s(PM-p). (9) 


Then we can use a and £ to describe the different mag- 
netic phases. The paramagnetic, ferromagnetic, antifer- 
romagnetic and ferrimagnetic phases correspond to the 
cases of a = 6 = 0, |a| > |8|, la] = 0,6 # 0 and 
0 Æ |a| < |8|, respectively. In the end of this paper, 
we will discuss the ferrimagnetic phase briefly. 

In terms of variables a and 8, the equations of motion 
for p and p° can be rewritten as 
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As in Ref. [22], our attention mainly focuses on the case 
near the critical temperature and we will ignore the dy- 
namics of p“ and p) since they are decoupled from 
equations of a and 8. To solve the equations of motion 
(10), we have to impose suitable boundary conditions 
at the horizon and AdS boundary, respectively. The two 
functions a and 8 should be regular at the horizon, which 
implies 
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at the horizon. So once the initial values of a and 8 are 
given at the horizon, one can integrate the equations (10) 
to get out the solution. The asymptotic behavior of the 
solution should not change the asymptotically AdS be- 
havior of the background geometry. As a result, near the 
AdS boundary we can neglect the nonlinear terms in (10) 
since the condensations for œ and 8 are small near the 
critical temperature, and obtain 
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The behavior of the solutions of equations depends on 
the value of m? + k +4. When m? + k +4 = 0, the 
asymptotic solutions will have a logarithmic term, we 
will not consider this case in the present paper. When 
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m? + k +4 Æ 0, we have the solution: 
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where a4 and 84 are all finite constants. To satisfy the 
Breitenlohner-Freedman (BF) bound in AdS4 space, the 
two model parameters m and k have to obey: 4|k| < 
17 + 4m?. On the other hand, as discussed in Ref. [22], 
in order to make the probe limit and the on shell action 
well defined, we should take 


6, >1, ô > 1, a, = b4 =0. (14) 


This can be understood in another way: when B = 0, the 
constants œ} and 84} should be viewed as the sources of 
the corresponding operators in the boundary field the- 
ory, according to the AdS/CFT duality. In order the 
symmetry to be broken spontaneously, one has to turn 
off the source terms. In addition, in order to make the 
p and p® condense and to realize the antiferromag- 
netic phase in the case without external magnetic field 
when the temperature is low enough, the model param- 
eters should violet the BF bound of 8 on AdS2, but not 
for a, which leads to 


3+2(m?+k)>0, 3+2(m?-—k) <0. (15) 


Considering (14) and (15) together, we see that m? and 
k should satisfy 


3 
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B. On shell action and free energy 


To see the phase transition of the model, we should 
calculate the free energy of the dual field theory, which 
is given by the on shell Euclidean action of the gravity 
theory. Since we are working in the probe limit, we need 
only calculate the part of the two polarization fields in the 
dyonic black brane background. A direct computation 
gives: 
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where A is the determinate of the induced metric at the 
AdS boundary and n” is the unit normal to the boundary. 
With the asymptotic behavior of the solution at r => oo 
and the restriction (16), we find that the boundary term 


in (17) vanishes. As a result, the free energy density for 
the part of p™ is 
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Here the total magnetic moment density is defined by 
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Here N, and N> can be interpreted as the magnetic mo- 
ments of two sublattices in the boundary. We can see that 
once the condensation happens, the nontrivial solution is 
always more stable than the trivial solution without con- 
densation, due to the fact that J < 0 and N > 0 as we 
will see shortly. 

The order parameter for antiferromagnetic phase tran- 
sition, i.e., the staggered magnetization, is defined as, 


N =N -N= -x | Ear (21) 


Because the free energy (18) needs to be invariant and 
magnetic field B will be changed into —B under the time 
reversal transformation, we can get the rules for magnetic 
moment and staggered magnetization under the time re- 
versal transformation, which is, 

p® Sol) Nis aN, (22) 
So if Ni Æ 0, then the time reversal symmetry is broken 
spontaneously. One should note that if the boundary is 
3+1 dimensional, magnetic components of polarization 
fields form two SO(3) axial vectors and staggered mag- 
netization is also SO(3) axial vector. Then if NÌ Æ 0, 
the time reversal and spatial rotational symmetries are 
broken spontaneously and simultaneously. This case will 
be further discussed in section IV. 


II. PHASE TRANSITION AND 
SUSCEPTIBILITY 


Now let us solve the equations (10) under the condi- 
tions (14) and (16). In what follows, we will take model 
parameters m? = —2,k = 1 as a typical example and 
work in grand canonical ensemble where the chemical 
potential is fixed. First, we will change the chemical po- 
tential to get different temperature of the black brane 
for any given magnetic field B. Second, we find out the 
suitable initial values of a and 8 at the horizon to sat- 
isfy the boundary conditions (14) through the shooting 


method. Third, note that the system has the following 
scaling symmetry: 


r— ar, (t,£, y) > a l(t, £,y), T > aT, u> ap, 
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Thus once obtain the solutions of (10), we can use the 
scaling symmetry (23) to find the corresponding solutions 
with the same chemical potential. 


A. Critical temperature and phase transition 


Let us first examine whether the polarization fields can 
condense in the antiferromagnetic manner when B = 
0. In order to do this, we will compute the values of 
(a2. + 87)~* for the different initial values ap and £o at 
the horizon and plot them in the plane of ao-89. Each 
singularity in this plane corresponds to a solution of (10) 
satisfying the condition (14). Because of symmetry, we 
only need to check the region of 69 > 0. The figure 3 
shows a typical example in the high and low temperature 
cases, respectively. In the high temperature case, we see 
that there is only one trivial solution locating at original 
point, which corresponds to a paramagnetic phase with 
a(r) = 6(r) = 0. When the temperature is low enough, 
we find there exist two nontrivial solutions which locate 
at Bo ~ 0.754 and ap œ 0, respectively. They lead to 
the solution with G(r) 4 0, a(r) = 0, which corresponds 
to an antiferromagnetic phase. 

Thus we see that the model indeed can give rise to a 
paramagnetism/antiferromagnetism phase transition in 
the case without external magnetic field. Since near the 
critical temperature, the value of (r) is very small and 
a(r) = 0, we can find the critical temperature by solving 
the linearized equation for (3: 
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Solving equation (24) with the condition of 6, = 0 
by shooting method, we find the critical temperature is 
Tn / œ 0.00925 for the case with chosen model param- 
eters. 

Now we can conclude that the system is in the para- 
magnetic phase with magnetic disordered when the tem- 
perature is higher than Ty and in the antiferromagnetic 
phase with two opposite magnetic moments when the 
temperature is lower than Ty. At Ty, there is a para- 
magnetism/antiferromagnetism phase transition. In ad- 
dition, since the magnetic moments of the two polariza- 
tion fields obtains expectation value, the time reversal 
symmetry is broken spontaneously. The numerical re- 
sults shows that this phase transition is a second order 
one with Ni = Nı — No œx \/1—T/Ty, where N, and 
Nə are the resulting magnetic moments of the two po- 
larization field condensations, respectively. As discussed 
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FIG. 3. The value of (af + 87)7' in the plane of ao-Go in 
different temperature. Top panel: T/u œ 0.0111. Bottom 
panel: T/u ~ 0.0083. 


in the end of subsection IIB, this corresponds to a time 
reversal symmetry breaking spontaneously. In addition, 
let us stress here that if the boundary spatial dimension 
is three, the spatial rotational symmetry is also broken 
spontaneously, since the nonvanishing magnetic moment 
chooses a direction as special. These results are consis- 
tent with the mean field theory description of the para- 
magnetism /antiferromagnetism phase transition given in 
appendix A. 


B. Static susceptibility in the paramagnetic and 
antiferromagnetic phases 


Two opposite oriented magnetic moments are one of 
the characteristic properties of antiferromagnetic mate- 
rial. Another remarkable one is the behavior of suscep- 
tibility density of the material in the external magnetic 
field. The static magnetic susceptibility density is defined 
by 

. ON 

X= m OB’ 2) 
When we turn on the external magnetic field B, the func- 
tions a and £ are both nonzero in any temperature. In or- 
der to compute the susceptibility density defined by (25), 
we need to shoot for the boundary conditions (14) with 
two parameters for equations (10) under the given exter- 
nal magnetic field B. This make us have to find the zero 
points of a nonlinear function with two variables. How- 
ever, since by definition, x involves only the behavior of 


B — 0, the problem can be simplified in the following 
way. 

In the case of T > Ty, taking the result a = 6 = 0 
when B = 0 into account, we expect that a(r) is in the 
same order as B, while 8 is in the order as B?. So in 
the case with weak magnetic field, we need not consider 
the equation for 8 and can neglect the nonlinear terms 
in the equation of a. In that case, we just need to solve 
the linear equation: 


tad 7 2 
n fla 2f 4 m*+k 2B 
=0. (26 
Vip (Teta + ar e ap 0 
Because a x B, we can set B = 1 simply. Solving equa- 


tion (26) with the boundary condition a, = 0, the sus- 
ceptibility density then can be computed by 


a 


x=N=-¥ | l ar. (27) 


In the case of T < Ty, the system is in the antifer- 
romagnetic phase. When B = 0, the solution of (10) 
has a = 0, 6 £0. When turn on a small B, we expect 
that a(r) is still in the same order as B, which leads to a 
term with the order of B? in the equation of 8. However, 
since the calculation of the susceptibility density only in- 
volves the linear term with B, thus equations (10) can 
be approximated reasonably to 


meg | JB? (74 4 +t *) 5-0, (28) 
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and 
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(29) 


In that case, the equation of 8 is decoupled from a. Thus 
we can solve 3 firstly and then put it into equation (29) 
and obtain a by solving (29). We can set B = 1 in 
(29) once again because a x B. Then the susceptibility 
density can be computed also through (27), once we have 
the solution a(r). 

Figure 4 shows the behavior of the susceptibility den- 
sity near the critical temperature Ty. In the paramag- 
netic phase, we can see that the susceptibility density 
increases when the temperature is lowered. Near the 
critical temperature, the susceptibility density satisfies 
the Curie-Weiss law of the antiferromagnetism (2). Con- 
cretely, for the chosen model parameters, we have 


à? / ux  0.0827(T/Tn + 13.65), (30) 


where 0/u ~ 13.657, /u = 0.1263. This is consistent 
with (2). When T < Ty, the susceptibility density in- 
creases when the temperature increases. The susceptibil- 
ity density has a clear peak at the critical temperature. 
These results are in agreement with the ones shown in 
figure 2 and with the results from the mean field theory 
discussed in appendix A qualitatively. 
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FIG. 4. The susceptibility density versus the temperature. 
(Top panel) The behavior of the inverse susceptibility den- 
sity in the paramagnetic phase near the critical temperature. 
(Bottom panel) The susceptibility density in the antiferro- 
magnetic phase and paramagnetic phase. 


C. Influence of strong external magnetic field 


As mentioned in section I, under the temperature be- 
low Néel temperature Ty, once the antiferromagnetic 
materials are put in an external magnetic field, the sys- 
tem will show a net magnetism because the two kinds of 
magnetic moments cannot offset each other. If the ex- 
ternal magnetic field is so strong such as B > Be, the 
magnetic moment structures will perform as paramag- 
netic materials where both two magnetic moments in two 
sublattices are the same and parallel to the direction of 
external magnetic field. According to the definition on a 
and £ in (9), this phenomenon corresponds to the case 
ofa #0, 6 £0 with the nonzero magnetic field B < Be 
and a #0, 8 =0 if B > Be in our model. 

In order to check these, let us see equation (10) again. 
For a finite magnetic field B, we can see that the equation 
for a is an inhomogeneous equation. So for any tempera- 
ture, there is always a nontrivial solution a(r) 4 0 which 
satisfies the condition of (14). From equation (10), we 
can read off the effective mass of a and ( near the horizon 
as 


2 
2 =m? tk- VE + 2r f’, 
3Ja? ey) 
m3 =m? —k-— Ar4 + 2rf'. 


From these relations, we can see that the increasing of a 


will increase the effective mass of 3, which will suppress 
the condensate of 8. If B is not very large, then a(r) is 
also not very large so that 6 can still condense. Then we 
can get solution of (10) 


a#0, B40, 


However, there is a critical magnetic field Be, at which 
the effective mass of 6 becomes large enough so that the 
condensation of 8 disappears. 

Of course, the critical magnetic field Be depends on 
temperature T. In the zero temperature case, the AdS2 
BF-bound of @ is 


when T < Ty, and0<B< Be. (32) 


3 
34+2(m? —k-— 770) >0. (33) 


Here ao is the initial value of a at the horizon. when the 
external magnetic field B gets increased, then ag will also 
increase. As the bound (33) holds, the condensation of 8 
disappears. For the finite temperature case, the critical 
magnetic field can be obtained by solving the following 
equations 
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(34) 


with the restrictions (14) and the initial condition 
6|r,=1 = 1. Figure 5 shows the behavior of critical mag- 
netic field in different temperature. We can see that Be 
will increase when temperature decreases. This agrees 
with the physical expectation. When the temperature is 
lowered, the interaction between two magnetic moments 
is stronger and stronger. In order to destroy this interac- 
tion, the external magnetic field needs larger and larger. 
Near the Néel temperature Ty, we can fit the behavior 
of B? with respect to 1 — T/Ty. The numerical results 
show a square root form 


B? ~ 7.9 x 10°T3(1 — T/Ty). (35) 


Because the staggered magnetization continuously de- 
creases to zero when external field increases to Be, there 
is a second order phase transition at B = Be. One can 
confirm this second order phase transition by comput- 
ing the free energy at the critical magnetic field and find 
that the free energy and its first order derivative are con- 
tinuous at the critical point. The critical behavior near 
B = B. has been investigated recently in [40]. 


IV. 4+1 DIMENSIONAL CASE 


In the previous section, our calculation is focused on 
the 3+1 bulk dimensional case, which corresponds to a 
2+1 dimensional boundary theory. In that case, only 
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FIG. 5. The critical magnetic field versus the temperature. 


the nonvanishing magnetic moment comes from the com- 
ponent M;,, such that the spatial rotational symmetry 
is in fact not broken even in the condensed phase, al- 


though the time reversal symmetry is broken sponta- 
neously. Now we generalize the model into the 4+1 
dimensional case in order to obtain a 3+1 dimensional 
boundary theory. In that case we can see clearly that 
the spatial rotational symmetry will be broken once the 
condensation happens. For this, we only need change the 
geometry of (6) into the form in the 4+1 dimensional 
case. In the 4+1 dimensional case, the condition (16) 
should be replaced by 


|[3+m?|<k<6+m?’. (36) 
In the 3+1 dimensional bulk case, the magnetic field is a 
pseudoscalar and the magnetic part in polarization field 
has only one component. The directions between these 
two components are parallel or antiparallel. In the 4+1 
dimensional bulk case, the components of magnetic parts 
of two polarization fields compose two SO(3) vectors, can 
this model also give two antiparallel magnetic moments 
rather than a non-collinear arrangement? The answer is 
yes. A brief proof is given as follows. 

In general, the ansatz for polarization fields can be 
written as 


Mi) =— p® (r)dt A dr + po) (r)da A dy + pe) (r)dz A dx 
+p) (r)dy Adz, a=1,2. 


(37) 


When the external magnetic field is absent, the equations 
of motion read 


a (b) 
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| 
where (a,b) = (1,2) or (2,1) and i = {a,y,z} and at the horizon remain, which are, 
(a)2 _ „(a)2 (a)2 (a)2 
pO? = pz + pz” + py’. Because of the global SO(3) lc 0)(rp), 02)(rx)}. (41) 


symmetry of the system, by adjusting the directions of 
axis {x,y,z}, we always can set following initial values at 
the horizon © 


PP (rn) = PP Cr) = OF (ra) = 0. 


By regular requirement at the horizon, the initial values 
of {p} have following relations 


(39) 


Ok + pm? 2, J pa? pi) (40) 
ArT : 


po = 2 + 


Combining the expressions (40) and (39), we have 
os (rp) = o (rn) = 0 and only three free parameters 


6 In fact, the global SO(3) symmetry permits us to set this condi- 
tion at any point. But one cannot set this condition in an open 
interval or two different points in general. 


Once given the values of them, we can integrate the equa- 
tions (38) to the boundary to obtain the whole solution. 
But near the boundary, by the variable substitutions sim- 
ilar to (9), the linearized equations give following asymp- 
totic solutions, 


1 1 
pi? = zi + Bi), pP = 3 (a — ĝi), 


Z Z 
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(42) 


As mentioned in the previous section, we need to impose 
following boundary conditions, 
Qi+ = i+ = 0, t=2,Y, 2, (43) 


which give 6 boundary constraints. However, we have 
only three free variables (41) at the horizon. As a result, 


in general, there is no solution except for the case that 


py 


components p (r), as discussed in the previous section, 
under the restrictions (36), they will condense with the 
same magnitude, but opposite signs. 

Thus we conclude that our model can indeed realize 
the antiparallel magnetic moment structure in the mag- 
netic ordered phase in the 4+1 dimensional bulk theory, 
which gives an expected 3+1 dimensional antiferromag- 
netic theory. Once the temperature is lower than a crit- 
ical value, the system will enter into antiferromagnetic 
state with time reversal and spatial rotational symme- 
tries broken spontaneously and simultaneously. 


(r) = ps?) (r) = 0. As to the remaining nontrivial 


V. SUMMARY 


In this paper, we have presented a holographic model 
to realize the paramagnetism /antiferromagnetism phase 
transition in a dyonic black brane background by intro- 
ducing two real antisymmetric polarization fields in the 
bulk. These two polarization fields couple to the back- 
ground gauge field strength and interact with each other. 
There exists a critical temperature Ty in this model if 
the model parameters are in some suitable region. 

In the case without the external magnetic field, these 
two polarization fields are both zero and the condensa- 
tion does not happen in the high temperature region of 
T > Ty, which corresponds to a paramagnetic phase. 
However, when the temperature is lowered to the region 
of T < Ty, the condensation happens spontaneously 
and the resulting two magnetic moments appear in an 
antiparallel manner with the same magnitudes, which 
leads to an antiferromagnetic phase, where the time 
reversal and spatial rotational symmetries are broken 
spontaneously. As a result, this is a kind of paramag- 
netism/antiferromagnetism phase transition at the criti- 
cal temperature T = Ty. 

When a weak external magnetic field is turned on, 
the resulting magnetic moment of the system is always 
nonzero. We discussed its response to the external mag- 
netic field by calculating the magnetic susceptibility den- 
sity. The results show that the magnetic susceptibil- 
ity density has a peak at the critical temperature and 
satisfies the Curie-Weiss law in the paramagnetic phase, 
which agree with the properties of susceptibility density 
in a realistic antiferromagnetism. 

Although most calculations are made in 3+1 bulk di- 
mensions, our results can also be extended into the 4+1 
bulk dimensional case in order to realize a 3+1 dimen- 
sional boundary theory. In details, in the 4+1 bulk di- 
mensional case, we proved our model can realize the time 
reversal and spatial rotational symmetries broken spon- 
taneously and simultaneously and the two magnetic mo- 
ment vectors are indeed antiparallel. 

This model can be generalized to describe the ferrimag- 
netism/paramagentim phase transition. Note that in the 
ferrimagnetic phase there exist two antiparallel magnetic 
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moments, but with different magnitudes. This can be 
archived by making the two polarization fields have dif- 
ferent masses, which lead to these two fields are not con- 
densed at the same temperature, or/and by considering 
different self-interacting constant J, which leads to the 
induced magnetic moments having different temperature- 
dependent behaviors. In such a model, we expect to have 
a complex spontaneous magnetization behavior and var- 
ious susceptibility curves. We are going to study these 
issues in the future. 

In appendix A, we give a mean field theory description 
of the antiferromagnetic phase transition, which can be 
regarded as a comparison. One may wonder whether it 
is worth building a holographic model to describe such 
a phenomenon, because it seemingly looks that what 
one can obtain from the holographic model can all be 
archived from a suitable mean field theory model. How- 
ever, essentially these two kinds of descriptions are differ- 
ent. First, the instability in the Landau-Ginzburg mean 
field theory must be put in by hand, while it arises nat- 
urally in the holographic setup. Second, the Landau- 
Ginzburg model is only valid near the transition point, 
whereas the gravitational description can characterize 
the whole dynamics. For a given bulk action, scanning 
through values of model parameters corresponds to scan- 
ning through many different dual field theories. In that 
sense, a simple holographic model is of some universal- 
ity, i.e., the results may be true for a large class of dual 
field theories, and are quite insensitive to the details of 
their dynamics. Third, in a top-down holographic setup, 
in principle one can know the details of dual field theory, 
while it is impossible in the mean field theory description, 
although our holographic model is a bottom-up one. 

There are many critical phenomenons involving strong 
correlation in condensed matter physics, which are con- 
sidered to be controlled by spontaneous magnetization 
or something relevant to it. The toy model presented 
in this paper can be regarded as a new starting point 
to investigate them from the holographic viewpoint. In 
particular, it gives us a wide possibility to investigate 
the coexistence, competition or interaction between mag- 
netic orders and superconducting orders, by combing this 
model with the holographic superconductor model [41]. 
The highly possible and very important applications are 
to investigate the coexistence and competition in ferro- 
magnetic superconductor or superconducting ferromag- 
netism [42—45] and the antiferromagnetic state in uncon- 
ventional superconductor [46]. All these works are ex- 
pected to be realized in the future. 
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Appendix A: Mean field theory for 
antiferromagnetism 


In this appendix we give a brief review on the mean 
field theory description for antiferromagnetism. Accord- 
ing to magnetic moment structure shown in figure 2, one 
can introduce two magnetic moments Mı and M2. Near 
the critical temperature, the free energy can be assumed 
to be of the form 


TM? H M2) + k Mı Mə cos 0 


(A1) 


b — > > > 
+34 + M$) - H- M, + Mə), 


where a, k and b are all positive constants, 0 is the 
included angle between two magnetic moments. is 
the external magnetic field, and Gop stands for the ir- 
relevant part. Comparing with the model of ferromag- 
netism/paramagnetism phase transition [28], one can 
view this free energy as the one of two magnetic moments 
with interaction between them. If set M 1 = Mo, we 
then get a ferromagnetism/paramagnetism phase transi- 
tion model with a negative critical temperature. Using 
the conditions of thermodynamic equilibrium, we have 


bM? +aT Mı + kMə cos 0 — H cos ¢, = 0, 
bM3 + aT Mə + kM, cos6é — H cos ¢2 = 0, 
Mı Mə sin = 0. 


(A2) 


A > 
Here ¢, and ¢2 are included angle of H and Mı and M2 
respectively. At first, let us consider the case of H = 0. 
The solutions for the equations (A2) are 


> > 
CER T > Ty = k/a, 


Mi =-Mə, Mil = VO- T/TN)k/b, T < Ty. 
(43) 
Thus we can see that when the temperature crosses Ty 
from the high temperature regime T > Ty, the system 
becomes magnetic ordered from disordered with a pair 
of opposite orientated magnetic moments. As the mag- 
netic moments condense in an opposite parallel manner, 
the material does not show any macroscopic magnetic 
moment. 
We can compute the difference of free energy between 
the antiferromagnetic phase and paramagnetic phase 
when T < Ty. A direct calculation gives 


T < Ty. (A4) 
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So when T < Ty, the antiferromagnetic phase is more 
stable than the paramagnetic phase. 

When the external magnetic field appears, using equa- 
tions (A2) and solutions (A3), we can obtain the mag- 
netic susceptibility for the case with the two magnetic 
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FIG. 6. (Top panel) The susceptibility versus temperature. 
We can see a peak at the critical temperature T = Ty. (Bot- 
tom panel) The free energy of the model in the high temper- 
ature regime T > Tyn and low temperature regime T < Ty. 


moments being (anti-)parallel to the magnetic field as 


2 
_ gp (2M _ 3M, _ ) KON +1)’ E EN 
XI = AoH t am) T 1 ee 
ti 0 
(A5) 


In figure 6, we plot the susceptibility (A5) and the free 
energy of the model in the case with k = a = 1, b = 1/3 
and Go = 0. One can find that the susceptibility satis- 
fies the Curie-Weiss law (2) in high temperature regime 
T > Ty. At the critical temperature, the susceptibility 
has a peak. When T < Ty, the susceptibility decreases 
when the temperature lowers. On the other hand, the 
free energy shows a spontaneous breaking of time rever- 
sal and spatial rotational symmetries when T < Ty. As 
a result the model of (A1) describes a second order phase 
transition and the order parameter is the staggered mag- 


> > > 
netization Mt = Mı — Mo. 
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